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Subroutine  ISOPLT  was  developed  for  the  computer  produced  di-play 
of  a  solution  space  surface.  The  solution  space,  for  the  sample 
plots  given  in  Appendix  B(  is  that  of  various  flow  functions  of  a  time 
dependent,  viscous,  incompressible  fuild  flow  for  various  boundary 
configurations.  Also  shown  in  Appendix  B  is  the  surface  plot  of 
SIN  (8* (X- 1 ) /XL+ 1/ 4 ( Y-l ) } +1 . 0  ^ 

Appendix  A  gives  a  sample  of  the  calling  sequence  for  equation  (1) 
in  the  MAIN  program,  and  also  a  listing  of  subrouting  ISOPLT. 

The  algorithm  used  in  ISOPLT  can  be  stated  as: 

1.  Draw  the  right  most  line  first,  moving  from  right  to  left. 

2.  Lift  the  pen  (blank  the  display  vector)  when  a  line  segment 

drops  below  (vertical  reference)  any  previously  drawn  line  segment, 
when  the  line  segment  again  moves  above  the  vertical  reference  line.) 

PROGRAM  VARIABLES 

TA ( I , J )  -  functional  value  at  the  point  I,J 
M  -  number  of  columns  in  TA 
N  -  number  of  rows  in  TA 

XP  -  size  of  plot  (in  inches)  in  horizontal  direction 

YP  -  size  of  plot  in  vertical  direction 

BND(I,J)  -  mode  point  type 

BND  ( I , J )  =  0  interior  point 

=  1  boundary  point 

=  2  exterior  point 

A ( I , J )  -  spacing  between  nodes  at  (I.J) 
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A  ( 1+  i ,  J  )  =  distance  from  A(I,J)  to  A(I+1,J) 
A  ( I , J  +  l )  =  distance  from  A(I,J)  to  A(I,J+1) 


etc . 

In  the  program  for  which  ISOPLT  was  developed,  both  A(I,J)  and 
BND  f I , J )  were  calculated  from  the  intersection  of  a  physical  boundary 
description  and  a  mesh  grid.  In  general,  however,  the  arrays  A(I,J) 
and  BND ( I , J )  may  be  initialized  as  shown  in  the  sample  MAIN  program 
in  Appendix  A. 
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(3 L  FOR  MIN 

DIMENSION  TA(61,26) ,  BND(61,26),  A(61,26) 
INTEGER  END 

COMMON  /D1  /  CAL. CMP  ,  TYPE 
DATA  CAL CMP/ 'CAL CMP’/ 

DATA  iype/'calcmp'/ 

M  -  61 

N  =  26 

CALL  I DP LOT 

CALL  PLOT  (.3  O', 6  .0,-3) 

DO  1  L  =  1  ,M 
DO  1  J-l  ,N 

X  -  8 . 0*FLOAT(1-1) / FLOAT (M-l) 

TA(I,J)  »  SLN (X  t  0,25* FLOAT ( J -1 ) )  +  1. 
AU.J)  ^  1-0 

1  BND(1,J)  =  0 

DO  2  L=1 ,M 
BND( 1 , J )  =  1 

2  BNDU.N)  =  1 

DO  3  J-l ,N 
BND(l.J)  =  1 

3  BND(M.J)  =  1 

XP  o  0 

yp  -  2,5 

CALL  ISOPLT (M,N ,XP ,YP ,TA,BND,A) 

CALL  PLOT (9  - 0 ,-3 .0 ,-3) 

CALL  F1N1 
CALL  EXLT 
END 
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@1  Fu  SOP! T 

SUBROUTINE  ISOPLTIM.N ,XP ,YP ,TA ,BND ,A) 

Di MENS  ION  BND  (61,26),  A(61,26) 

DIMENSION  TA(6 1 , 26)  ,  U(75),  V(75),  T(75),  S(75)  ,  TEM(75) 
DIMENSION  UA( 75) ,  VA( 75) 

REAL  Ml,  M2 
LOGICAL  DOWN,  FIRST 
DIMENSION  FI,(  3000) 

COMMON/ PL  1/  SCALE,  Ft 
COMMON  /DI /  CAL CMP ,  TYPE 
IN  TEGER  BND,  FL, ,  TYPE,  CALCMP 


N  t  -  N  *■  ) 

N2  -  N  +  2 

DY  80*YP/ FLOAT (N-l) 

DX  =  DY 
X  =  XP  +  DX 
ANG  -  0  8bb02540 
ER  =  0  2  *ANG*DX 
ZSCALE  -  0  0 
DO  1  JM  ,N 
DO  1  l  —  1 ,  M 

LF ( BND( I , J ) , EQ  7)  GO  TO  1 
ZSCAI.E  AMAX)  (ZSCALE,  ABS  (TA(  I ,  J)  )  ) 
1  CONTINUE 

DO  10  II  -  1  ,M 
J  -M+l-IL 

X  =  X  -  DX 

Y  =•  'DY 
0(1)  -  ANG*X 
V(  1)  --  -  5*X 
UA(1)  -0.0 
V'A(l)  *  0-0 
JJ  -1 

NN  -  N  -  I 
DO  2  J  =  i  ,NN 
J  J  -  J  J  +  L 
UA(JJ)  -  0.0 
VA(JJ)  =00 

Y  =  Y  +  DY 

U(JJ)  =  ANG  MX  +  Y) 

LF ( BND( I ,J ) , EQ • 0)  GO  TO  31 
I F ( BND( 1 , J ) .EQ  1)  GO  TO  30 
V(J.J)  --  5MY-X) 

I F ( J  EQ  N)  GO  TO  2 
I  FlBNDd  ,J+)  )  EQ.2)  GO  TO  2 
UA(jJ)  -  ANG*DY* (l.-A(l,J+l)) 

VA(JJ)  5*DYM  1 .  -A(  I  ,J  +  1)) 

J  j  J  J  +  1 
U(J.J)  -  U(JJ-l) 

V'.JJ)  =  .  5* ( Y-X)  +  TA(1,J  +  1) /ZSCALE 
UA(jJ)  =  UA(JJ-l) 

VA(JJ)  =  VA(JJ-l) 

GO  TO  2 

31  V(JJ)  =  5* (Y-X)  +  TA( I ,J+1 ) /ZSCALE 
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UAUJ)  ANG*DY*U  .-AU,J+D) 

jj)  5*DY*U  -A(I,J  +  1)) 

IHj  EQ.NN)  GO  10  2  2 
IFlBNDU  ,Jt)  )  EQ  0)  GO  TO  2 
1 F (.  J+  2  GI  i N )  GO  10  2 
l Ftfv  NDt 1 ,Jt 2 )  LT ,2)  GO  TO  2 

22  J  I  -  JJ  t  1 
L'tJJ)  =  U(JJ-l) 

V  C  J  J  )  =  5*(  Y-X) 

UAlJJ)  =  UA(Jj-I) 

VA(JJ)  '  .’aUJ-D 
GO  10  2 

30  LF(BNDU.Jtl)  EQ  2)  IA(I,J  +  1)  =  0,0 
V(JJ)  -  3HY-X)  f  TA(I,J  +  1)/ZSCALE 
UAUJ)  =  ANG*DY*(l.-A(ItJ+U) 

VAlJ  J  )  =  5*DY*(.1--AU,J  +  1)) 

1  F  ( J  I.T  ,NN  )  GO  TO  2 
Jj  =  JJ  +  i 
UUJ(  =  UUJ-i) 
vuj)  =  5*a-x> 

2  CONIINUE 
DOWN  =■  i  RUE 

1  F  (  J I  G'f  .))  GO  TO  19 

C  PLOT  FIRST  LINE  WITHOUT  HIDDEN  LINE  ALGORITHM.. 

CALL  PLOT  HU  ( 1 )  ,  V(l) ,  3) 

DO  21  J  =  2  ,JJ 

21  CALL  PLOT  1 (U ( J ) -UA( J ) ,  V(J)-VA(J),  2) 

DO  2  3  K  =  2  ,jJ 

23  T  (K)  =  -10  0 
GO  TO  33 

19  CALL  PLOIl(UU)  ,  V(l)  ,  3) 

DO  20  V--2  ,J 

20  CALL  PLOrHU(K)-UALK)  ,  V(K)-VA(K),  2) 

S(l)  1  Ul  1  > 

T(l)  »  Vll) 

FIRST  -  FALSE, 

i F( V(4 )  ,  L  T  TO))  FIRST  =  .TRUE. 

DO  8  K=4 ,JJ 
DO  4  3  KK=  1  ,N1 

IF(ABS(SLKK)-U(K))  GT ,ER)  GO  TO  43 
IFLV(K) -1 (KK) )  50,  12,  12 
50  I F ( DOWN )  GO  TO  14 

CALL  PLOT  HULK)  ,  V(K)  ,  3) 

DOWN  *  FALSE 
GO  TO  8 
4  3  CONTINUE 

PRINT  44,  K 

44  FORMAT (1H0,  30X  ,  8HHELP . . ,  ,  2HL=I3) 

RETURN 

12  FIRST  =  FALSE. 

13  I F (  NOT  DOWN)  GO  TO  14 

CALL  P LOT l (U(K)-UA(K) ,  V(K)-VA(K),  2) 

GO  10  8 

14  Ml  =  a  IKK)  -  T (KK-1 ) ) / (S (KK)  -  S(KK-l)) 
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IHABSlU(Kj-Ui'K-l))  .GT  .ER)  GO  10  4 

:>S  U(K; 
bb  10  5 

4  M2  =  <V(K)  -  V  (K-l) )  /  (U(K)  -  U(K-l)) 

So  (M2*U (K-l)  -  Ml *S (KK-1  +  I(KK-l)  -  V(K-1))/(M2  -  Ml) 

5  IT  ~  Ml  *  (SS  -  S  (.KK-1) )  +  T(KK-l) 

If  (at.  <K - 1 )  -ER  GT  .SS  .OR.  SS  .GT  ,S(KK)+ER)  GO  TO  16 

IE (DOWN)  GO  TO  27 
GALL  FL01HSS,  TT,  3) 

CALL  PLOT  1 (U (K) -UA(K) ,  V(K)-VA(K),  2) 

DOWN  -  TRUE . 

GO  TO  8 

2  7  CALL.  PLOT  KSS-UA(K)  ,  TT,  2) 

lb  CALL  PLOI 1  (U(K)  -l)A(K) ,  V(K)-VA(K),  3) 

DOWN  =•  FALSE. 

8  CONTINUE 

33  S ( 2  )  -  U(l) 

DO  40  K= 2  ,N 
DO  4  5  J  K,JJ 

IF'(UU)  GT . S (K)+ER)  GO  TO  40 
45  CONTINUE 

40  S(K+J)  =  U(J) 

C  STORE  THE  MAXIMUM  VALUE  OF  THE  V(K)  LINE  IN  ARRAY  T(K) 

T  ( 1 )  =  V(l) 

DO  41  K=2  ,N1 
IEM(K)  --  -10  0 
DO  4 2  KK=K , JJ 
IF (U ( KK) -S  CK) )  42,  51,  42 
51  TEM(K)  =  AMAXKV(KK),  T(K-l)  ,  TEM(K) ) 

42  CONTINUE 

41  CONTINUE 

DO  11  K-2  ,N 1 
11  T(K)  =  TEM(K) 

6  CONTINUE 
10  CONTINUE 

1 F ( TYPE , NE  CALCMP )  CALL  SNDFLE(FL) 

RETURN 

SUBROUTINE  PLOT  1 (A ,B ,KK) 

DIMENSION  FL ( 3000) 

COMMON  /D1 /  CALCMP,  TYPE 
COMMON /PLT7  SCALE,  FL 
COMMON /PLT2/  ISW,  1X1,  IY1 
INTEGER  TYPE,  CALCMP,  FL 

IF (TYPE  NE. CALCMP)  GO  TO  1 

CALL  PLOI (A ,B ,KK) 

GO  TO  6 

1  1 F (KK . LT  0)  GO  TO  6 
IX  =  /O.O’-A  +  300.0 
IY  =  70  0*B  =  600.0 
GO  TO  (2  ,2,3)  ,  KK 

2  IF  (KK  .  EQ  ISW.)  GO  TO  7 
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(".ALL  INHPUO) 

i  !  1NE(H.  ,  1 A  ) 

CALL  INnPLlV) 

I5W  - 

(.All  I  ;  NE(  L  L  ,  IX  ,  1  V  ,  I  A) 
(jO  I U  A 


APPENDIX  B 


